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We discuss the effect of elastic deformations on the electronic properties of bilayer graphene 
membranes. Distortions of the lattice translate into fictitious gauge fields in the electronic Dirac 
Hamiltonian which are explicitly derived here for arbitrary elastic deformations. We include gauge 
fields associated to intra- as well as inter-layer hopping terms and discuss their effects in different 
contexts. As a first application, we use the gauge fields in order to study the recently predicted strain- 
induced Lifshitz transition for the Fermi surface at low energy. As a second application, we discuss 
the electron-phonon coupling induced by the fictitious gauge fields and analyse its contribution to 
the electrical resistivity of suspended bilayer membranes. Of special interest is the appearance of a 
linear coupling for flexural modes, in stark contrast to the case of monolayer graphene. This new 
coupling channel is shown to dominate the temperature-dependent resistivity in suspended samples 
with low tension. 
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I. INTRODUCTION 

Graphene, a monolayer of carbon atoms, is the 
only two-dimensional (2D) conducting elastic membrane, 
where electrons have been shown to behave as mass- 
less Dirac fermions.^'^ Since the experimental realisa- 
tion of graphene, several authors investigated the in- 
terplay between its electronic properties and externally 
induced mechanical deformations (see, e.g. Ref. [3] and 
references therein). This question acquires a particu- 
lar interest since elastic distortions of the honeycomb 
lattice translate into fictitious gauge fields in the elec- 
tronic Dirac Hamiltonian.^ ^ Indeed, in presence of de- 
formations, the low-energy massless Dirac Hamiltonian 
of electrons in monolayer graphene is described by = 
va - {p -\- eA) with the fictitious vector potential A ex 
[u^^ — Uyy, — 2'^^,^),'^''^ expressed in terms of the strain 

tensor u,, of the 2D membrane. The latter reflects 

''J 

the symmetry of the system with respect to the plane 
of graphene, as it is linear in the in-plane deforma- 
tions and quadratic in the out-of-plane (flexural) distor- 
tions. As a direct consequence of this vector potential, it 
has recently been observed experimentally that strained 
graphene bubbles exhibit the formation of Landau levels 
corresponding to magnetic fields as large as 300 T.^ This 
adds significant motivation to recent works devoted to 
the investigation of the so-called strain-engineering, i.e. 
the tailoring of the electronic properties of graphene via 
controllable elastic deformations 

Bilayer graphene membranes show an even richer sce- 
nario which is attracting ever growing interest. This is 
partly due to their larger potential for device applica- 
tions, since a gap in the spectrum can be induced by 
external electric fields .^^^^^ From a fundamental point 
of view, the low-energy electronic band structure of bi- 
layer graphene shows fascinating structures, such as the 



appearance of four massless Dirac cones which evolve 
into a massive quasiparticle spectrum at higher energy.El 
Given the remarkable tunability of the electron density in 
graphene devices, this opens the possibility of observing 
a phase transition in the topology of the Fermi surface 
(FS) as a function of the electron doping, the Lifshitz 
transition (LT).^^ 

So far, the investigation of the effects of elastic defor- 
mations in bi layer graphene has been limited to uniaxial 
defor mat ions .^^^^ In particular, the form of the fictitious 
gauge fields associated to arbitrary distortions is still un- 
known. These fields can result in dramatic effects on the 
electronic band structure, and are the starting point for 
any systematic investigation of electron-phonon coupling 
and of strain engineering in bilayer samples. In this paper 
we address this open question by deriving the form of the 
gauge fields and discussing how they affect the electronic 
Dirac Hamiltonian. 

As a first application of our analysis we consider the 
special case of elastic deformations yielding uniform fic- 
titious gauge fields (like uniaxial strain, pure shear, or a 
rigid shift between the two layers), including some that 
have been recently discussed in the literature.'^^^ These 
deformations dramatically modify the electronic bands 
at low energy and induce the annihilation of two mass- 
less cones beyond a critical value of strain. This affects 
the nature of the LT and the energy at which it occurs, 
with observable consequences on the formation of Landau 
levels (LL) under perpendicular magnetic fields and on 
the single particle density of states (DOS) at the Fermi 
level. The results of this section reproduce those that 
appeared recentl}/^^ and serve as a test for the general 
gauge fields we deduce. 

As a further application of the fictitious gauge fields, 
we discuss the contribution to the electronic resistiv- 
ity of suspended bilayer membranes due to electron- 
phonon coupling. In contrast with the case of mono- 
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layers, the two layers are not equivalent and the system 
does not exhibit a symmetry with respect to out-of-plane 
deformations. This is reflected in the appearance of a 
linear intrinsic coupling for flexural phonons in bilayer 
graphene, in contrast to the purely quadratic one present 
in monolayers .I^I^^H^ This linear coupling channel results 
in the dominant electron-phonon contribution to the re- 
sistivity for suspended samples with low tension. Our in- 
vestigation provides the basis for the microscopic study 
of the electro-mechanical properties of bilayer graphene 
membranes, including the fictitious magnetic fields re- 
sulting from arbitrary elastic mechanical deformations. 

The structure of the paper is as follows: In Sec. [ll| 
we briefly summarise the electronic properties of bilayer 
graphene, considering all hopping parameters in the orig- 



inal tight-binding Hamiltonian. In Sec. Ill we derive the 
explicit expression for the fictitious gauge fields in bi- 
layer graphene under arbitrary elastic deformations, and 
deduce the effective low-energy Hamiltonian for the two 
bands close to zero energy. Section IV summarises the 



electronic band structure in the absence of deformations. 
In Sec. [V| we introduce elastic deformations resulting in 
uniform fictitious gauge fields and discuss their effect on 
the electronic low-energy spectrum. We analyse the dis- 
placement of the massless Dirac cones and their anni- 
hilation under a critical strain. As a consequence, we 
explore the tailoring of the band-structure via control- 
lable elastic deformations, finally resulting in the possi- 
bility of inducing topological LT at the Fermi level. In 
Sec.|Vl]we discuss the electron-phonon coupling induced 
by the fictitious gauge fields in suspended bilayer mem- 
branes and analyse its contribution to the temperature- 
dependent component of the resistivity. Finally, in Sec. 




FIG. 1. Atomic structure of perfect bilayer graphene in Bernal 
stacking (top view) . The honeycomb lattice with full (dashed) 
lines corresponds to the upper (lower) layer / = 2 (/ = 1). The 
atomic sites A2 and Bl coincide once projected on the plane. 
The three vectors e^, {j = 1,2,3) connecting A and B sites 
are indicated (see text for more details). 



tonian in the real space 
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VII we present our conclusions. 



^B2 i = l 



(1) 



II. ELECTRONS IN BILAYER GRAPHENE 

The electronic properties of ideal bilayer graph ene have 
been studied in the past by several authors. The car- 
bon atoms in the two layers are identified by a layer 
index (/ = 1,2, for the first and second layer, respec- 
tively) and a sublattice index {s = A, 5), to distin- 
guish inequivalent atomic sites. In each layer, carbon 
atoms form honeycomb lattices with nearest neighbour 
distance a = 1.42 A. The three nearest neighbours of 
type A are obtained from a central B atom by the dis- 
placement vectors ei = a (0, — 1), e2 = a (a/3/2, 1/2) and 
ea = a (—a/3/2, 1/2). In the Bernal stacking configura- 
tion for a bilayer in the x-y plane, the atoms in layer 2 
and sublattice A (i.e. of type A2) are located directly 
above those of type Bl, at a distance c c:^ 3.34 A, as in- 
dicated in Fig. [l] 

On the contrary, B2 atoms are not located above Al 
ones. In fact, the pairs A1-A2, B1-B2 and A1-B2 are 
all separated by the distance c = Vc^ + — 3.63 A. 
The starting point for the investigation of the electronic 
properties of bilayer graphene is the tight-binding Hamil- 
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-^A1,A2 Y Y I^A2)(R'A2 + - CZ| 

^A2 i = l 
3 

-^B1,B2 Y Y I^B2)(R'B2 + - CZ| + h.C. 

^B2 i = l 

where R^^ is the position of an atom in sublattice s and 
layer /, and |Rg^) is the ket associated to the correspond- 
ing localised orbital. For the interlayer hopping terms we 
used the identities R^2 = ^ 

^Al — ^A2 + ~ ^Bl — ^B2 

the tight-binding Hamiltonian, t^i ^/^/ represents the hop- 
ping energy between the two neighboring sites at R^^ and 
R^/^/. In bilayer graphene samples these hopping terms 
are given by t^j j^j = tj 2.47 eV (with j = 1,2)^ 

0.39eV,'^t^i^^2 = 73 - 0.315 eV,E3] and 
74 0.044 eV.^^ For the sake of gen- 
erality, in the Hamiltonian above we consider the case 
in which the intralayer hopping energies and ^2 can 



■A2 ~ -^^Bl^^^' -^B2 ~ R^i+e^ +CZ, 

e^- — cz. In 



^A2,B1 — 
^Al,A2 — ^B1,B2 
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be diflFerent. This can be of relevance for bilayers on a 
substrate, where the direct contact with the latter may 
affect the hopping energy in one layer with respect to 
the other. Electrons in the bilayer lattice are described 
by Bloch states of the form 



with the amphtude of the wavefunction on the sub- 
lattice s and layer I at wavevector k. In the 4x4 space 
of the Bloch amplitudes {u'^^\ u'^^\ u^^'' ) the 

Hamiltonian thus takes the form 



-74 /k 



-73/k -74/k -*l/k 

-t2/k -74/k 



-74/k 





-7 



-7 





(3) 



with = Xlj=iexp[zk • e^], resulting in four energy 
bands. Two of them are at high energy of order ±7, while 
the other two touch close to zero energy, in the vicinity 
of two Dirac points in the first Brillouin zone given by 
K^, with = r(47r/3V3a,0), where /^^ = 0. In the 
vicinity of the two Dirac points, the matrix Hamiltonian 
in Eq. (|3| is expanded as 
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(-) 
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\ -vip 
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V2P^ 



V4P' 
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-v^p 


-V2P 



vip^ \ 
v^p 

-7 



-V4P 

-V2P^ 



-7 



(4) 



-Vip \ 
-v^p^ 

-7 
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where p = p^ -\- ipy is the complex representation of the 
two dimensional momentum relative to the Dirac point, 
Vj = 3atj/2h for j = 1, 2 and Vj = 3a'yj/2h for j = 3,4. 

The two valleys are thus related by the symmetry i^o~^ = 



(+)*l 



We now introduce elastic deformations in 



the bilayer lattice and derive the consequent fictitious 
gauge fields in the Dirac Hamiltonian formalism. 



III. FICTITIOUS GAUGE FIELDS IN BILAYER 
GRAPHENE 



A generic elastic deformation in the bilayer graphene 
membrane induces a displacement of the atomic posi- 
tions which results in the modification of bond lengths 
between neighboring atoms. The corresponding change 
in the hopping energies results in corrections to the ma- 
trix elements of the 4x4 Hamiltonian Hq'\ yielding a 
shift of electronic momenta analogous to that stemming 



from a vector potential. As a result, mechanical defor- 
mations translate into fictitious gauge fields in the Dirac 
Hamiltonian. 

In addition, elastic deformations involving local varia- 
tions of areas induce fluctuations in the electronic den- 
sity that translate into scalar deformation potentials in 
each layer. ^ ^ As a result, symmetric deformations be- 
tween the two layers yield a global scalar potential in the 
bilayer Hamiltonian. In contrast, antisymmetric distor- 
tions yield a difference in the potential between the two 
layers that mimics the effect of an inter layer electric field. 
The latter mechanism thus results in a deformation- 
induced gap in the electronic spectrum even in the ab- 
sence of external gates. 

In monolayer graphene the fictitious vector poten- 
tial induced by mechanical deformations has been in- 
vestigated in the past.^ ^ The massless Dirac Hamil- 
tonian describing the low energy properties of elec- 
trons in a deformed suspended membrane has the form 
= va • (p + eA'^+)) in the valley r = +, 
with V the electron velocity, a = {a^^ay^a^) the 
pseudospin vector of Pauli matrices in the sublattice 
space and p = {p^^Py^O) the electronic momentum. 
The fictitious vector potential is expressed as A^+^ = 
ctq h/t (dt/da) {{u^^ — Uyy)/2^ ~'^xy) terms of the in- 
plane hopping energy t and of the strain tensor of the 2D 
membrane u^j = {d^Uj -\- djU^ -\- d^h djh) with ctq the 
unit matrix in the sublattice space. Here u(r) is the vec- 
tor field describing in-plane deformations and h{r) the 
scalar field associated with out-of-plane (flexural) dis- 
tortions. The symmetry with respect to the plane of 
graphene is reflected in the quadratic contribution from 
flexural deformations, as the effects on electrons due to 
the displacements /i(r) and — /i(r) are identical. The 
Hamiltonian in the other valley (r = — ) is given by 
= ^_ , leading to fictitious vector poten- 

tials with opposite signs in the two valleys. In particular, 
the fictitious magnetic fields associated to the vector po- 
tentials have opposite signs in the vicinity of the two in- 
equivalent Dirac points, since the effects of deformations 
do not break time-reversal invariance. 



A. Effects of deformations in the electronic 
Hamiltonian 

In order to calculate the vector potential in bilayer 
graphene, we analyse the shift of the atomic positions un- 
der a generic distortion. In-plane deformations in layer / 
are described by the two-dimensional vector field u*^^^(r), 
while out-of-plane distortions are associated to a scalar 
field /i^^^(r), so that the atom at position r is shifted 
to r + u*^^^(r) + zh''^\r). Thus, in the tight-binding 
Hamiltonian, the hopping term between |Rg^) and iRg/^/) 



t 



sLs'l' 



undergoes the change t^^ ^/^/ 



with 51^ 



Here 



Ls'l' 



u('')(R,,,,)-u«(R,,) + z[/.('')(R^,,,) 
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R 



- R 7 the variation of the correspond- 



ing bond length. Among the derivatives dt^i s'V l^^si , 



only the intralayer one dtj^^ / d£ 



'^ALBl 



-3t 



ALBl/^ALBl 



is known. For all the other cases we will assume typi- 



cal values dt^i^^^i^/di^i^^^i^ -Vsi,s'u'^si,s'U /^si,s'U ^ with 
r]gi of order one. Performing the Bloch band anal- 
ysis of the deformation-induced corrections to the tight 
binding problem, we can thus obtain the contribution to 
the Hamiltonian associated with elastic distortions. For 
this purpose it is convenient to introduce symmetric {S) 
and antisymmetric (A) deformations in the layer index 
as follows 



AS /A) 



(r) 



u(2)(r)±u(i)(r) 



(5) 



where the upper /lower sign is associated to the symmet- 
ric/antisymmetric channel. In terms of these modes, and 
using a long-wavelength expansion of the variation of the 
bond lengths J^^^^/^/, we deduce the corrections to the 
Hamiltonian in the two valleys (r = ±) in the Dirac ma- 
trix formalism as 
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(6) 



with 



^1=1,2 - 



gTi[u, 



(Oi 

ij J 



da 



,(0 



,(0 



(r) 



F — —2 
^~ dc 



2c dc 

2c dc 
dj 



(7) 



u 



(A)2 



Here li-P = {d^u^p 



rdjuf^ + d^h^^^ djh^^^)/2 is the strain 
tensor of the two-dimensional membrane in layer / and, 
in lowest order in the deformation fields, we find 

^[u(^), u(^), /i(^), /i(^)] = ac - iTd^h^'^) 

+2a(4^) -ir4^)) . (8) 



Thus the correction terms in the two valleys r = ± are 
related by the symmetry 6H^~'> = In Eq. ^ we 

also introduced the term Di representing the deforma- 
tion potential for the layer / associated to local variation 
of areas in a distorted elastic medium,^ with g the de- 
formation potential coupling constant. While the bare 



value of g is estimated around 20 — 30 eV,^ in graphene 
samples the coupling constant is reduced by screening so 
that it effectively depends on the electron density.l^^^ In 
contrast, gauge fields are not affected by screening.'^ As 
a result, one expects that the vector potential dominates 
in (significantly strongly) doped graphene, while the de- 
formation potential would be dominant in the immediate 
vicinity of the Dirac point. 

The terms Fi^\ ^ Eq. ^ are the same as those ap- 
pearing in monolayer graphene. They are linear in the 
in-plane deformations causing variations in the bond- 
lengths. However, the symmetry with respect to the 
x — y plane forces the coupling with out-of-plane deforma- 

(t) 

tions to be quadratic. In contrast, in the terms F3 ^ and 

^4^^ it is interesting to notice the appearance of a linear 
coupling between electrons and symmetric fiexural defor- 
mations {h'^^^) for the skewed interlayer hopping terms 
(A1-A2, 51-52 and Al-52). In view of the results for 
the monolayer, protected by the symmetry with respect 
to the plane, it may be surprising that such a term exists. 
However, due to the rotation between the layers involved 
in the Bernal stacking, the two layers are not equivalent. 
Thus, for skewed interlayer bonds, the symmetric fiexu- 
ral deformations h'^^^ and —h'^^^ induce different effects, 
resulting in the appearance of a linear residual contribu- 
tion in the fictitious gauge fields. In more detail, this can 
be understood by considering the bond lengths involved 
in the interlayer hopping terms different from the "verti- 
cal" Bl — A2 one. In Fig. [2] the interlayer skewed bonds 
are illustrated under a generic h^^"^ (and —h^^"^) defor- 
mation. The two are evidently different, due to the shift 
between the projected positions of the atoms involved. 
The analysis of the figure also reveals why the Fg^^ and 
^4^^ terms differ by the replacements h^^"^ —h^^"^ and 
u^^^ —u^^\ as these preserve the bond lengths in- 
volved in the corresponding hopping terms. 

As far as long wavelength antisymmetric fiexural de- 
formations h'^^^ are concerned, they correspond to a local 
modulation of the interlayer distance which preserves the 
structure of Bernal stacking. Their only effect is thus to 
induce a numerical renormalisation of the velocities 
and v^. In Eqs. ^ and ([8| this would merely lead to 
a sub-leading correction in h^^^ /c <C 1 and \p\a/Ti <C 1 
of the form 5tF = {2ac/?i) h^^^ ijVx + '^Py) 5 which can be 
neglected. 

Finally, the term F^ is associated to the vertical inter- 
layer hopping Bl—A2^ which does not involve any skewed 
bond. As a consequence it is not affected by symmetric 
deformations h^^"^ and u^*^^. The variation of the bond 
length stems uniquely from antisymmetric deformations. 
It is linear in the fiexural distortions h'^^^ and quadratic 
in the in-plane ones u^^^. 

The energy corrections fJ^^ (with j = 1,...,4) in 

6H^'^'> thus affect the unperturbed terms Vjp in Hq''^ 
as fictitious gauge fields acting on the electronic orbital 
degrees of freedom. The symmetry Hq ^ + 5H^~^ = 
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4' = 3" 



FIG. 2. Schematic description of the effect of deformations 
on bond lengths. This side view, taken along the direction 
Qj, shows the atoms involved in intra and interlayer hopping 
processes in a unit cell, a) Bilayer under flexural deforma- 
tions h^^^ and —h^^\ The difference between the hopping 
lengths 3' and 3^' (as well as between 4' and 4'') breaks the 
symmetry with respect to the plane and is responsible for the 
appearance of a linear coupling with h^^"^ in the gauge field 



Fg^^^ (and F]^^). The fact that 3' = 4" leads to the symmetry 
h^^^ -h^^^ between F^''^ and F^''\ b) same as in a), but 
for in-plane antisymmetric deformations u^"^-* and — u^^-^ . The 
equality 4' = 3" accounts for the symmetry u' 
between F^^^ and F^^K 



(A) 



(+)*| 



SH^^>, together with 



(-) 



veals that the fictitious gauge fields have opposite signs 
in the two valleys, as in the monolayer case. As a conse- 
quence, the fictitious magnetic fields generated by elastic 
deformations are also opposite in the two valleys, as re- 
quested by the fact that elastic deformations do not break 
time-reversal invariance. 



B. Explicit form of the fictitious vector potential 



In order to express the effect of deformations on the 
electronic momenta in terms of a vector potential, we 
can rewrite the total Hamiltonian H^^^ + 5H^^^ in the 
language of Pauli matrices acting on the layer space (E^, 
with a G {0, x, z}) and on the sublattice space {a^) as 



7- 



(9) 



H 



D 



2 



^0 ' 



with Sq and ctq the identity matrices in the corresponding 
spaces. The term collects all contributions linear 

in v-p and 5Hp'^^ the corrections F- (with j = 1,...,4), 
while collects the terms involving 7 and F^, and Hjj 
those related to the deformation potentials. In the latter 
we introduced the symmetric and antisymmetric compo- 
nents F>(^) = (F>i + D^)/2 and D^^^ = {D^ - D^)/2, 
corresponding to different variations of areas in the two 
layers. In analogy with the monolayer case, the Hamil- 



tonian Hp'^'^ can be written as 



/f(+)=vW-(p + eAW) 



(10) 



with Hq^J = V(+^ • p expressed in terms of the vector of 
velocity matrices 



v(+) 



with 

S., (g) a,.) + 



(11) 



2 -z^-y 



In parallel, for the term 5Hp^^ , by direct inspection one 
finds 

5Hl+)^eV^+)-A(+)=5H'^J^^+6H[il, with (12) 



5H, 



(+) 

Re,p 



Re[Fj 



(+)i 



Re[K 



(+)i 



(+) -Tv.rp(+)i5o±^ 



E„ (K) cr. 



Im[K 



(+)l ^0 - 



3 . 2 

leading to the vector potential 



- + Im[Fi+^] S„ ® (T, 



' 



(13) 
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V = 3at/2h this indeed coincides with the vector poten- 
tial of a monolayer. The coupling between the layers 
would still be present via the term which, however, 
does not involve a vector potential. 



In contrast to monolayers, for bilayer graphene in the 
general case (with v-^ '^3 7^ 0, 7^ and all the 

associated corrections ^j^^) the analytical expression for 
the vector potential is quite cumbersome, due to the ma- 
trix structure of the velocity vector V^+^. For illustra- 
tion purposes we consider the simpler case = '^2 — 



^ 1 — ^ 2 
compact expression 



and Va 



leading to the 



(14) 



The terms proportional to and ^3^^ are related to in- 



C. Effective low-energy Hamiltonian 

The 4x4 Hamiltonians H^^^ ^SH'^^^ in Eqs. (4| and ^ 
contain the complete information concerning the proper- 
ties of electrons and of their coupling to elastic defor- 
mations. The unperturbed Hamiltonian Hq'^ is diag- 
onalised in terms of four electronic bands, two of them 
touching at low energy, and two describing split modes at 
high-energies of order ±7, due to the "vertical" Bl — A2 
interlayer hopping. 

In order to analyse the low-energy sector of the spec- 
trum, we can produce an effective 2x2 Hamiltonian 
in the Al — B2 subspace along the same line as done, 
e.g. in Refs. 11 and |28l We consider the original 
4x4 Hamiltonian H^^^ = H^^^ + SH'^^^ as made of 
four 2x2 blocks H-j\ G 1,2) with the upper left 

(t) 

block H^-^^ describing the low energy sector. We intro- 
duce the matrix Green's function G^^"^ = (el — H^^^^ , 
with 1 the unit matrix, and by direct matrix inversion 



terlayer hopping processes and are associated with the we get Gy^ = ^el — Hy^ — {el — ^22^)~^^2i^) 



appearance of the Pauli matrices which involve 

mixing of the two layers. In contrast, the special case 
= Fg^^ = would lead to two decoupled layers as 

far as Hp'^^ is concerned, each characterised by the vec- 
tor potential eA^+) = ctq/v (Re[F^+)], Im[F^+)]). Us- 
ing the expression for F^~^^ in Eq. til and the velocity 



Thus the effective low energy Hamiltonian is obtained as 

= el - G[\^-' h[1^ - Hil\H^Y^Hil^ in the 
limit e <C 7. By keeping only the lowest non- vanishing 
order in Fj^"^ /j <C 1 and Djj <C 1, after lengthy but 
straightforward steps, we obtain the effective low-energy 
Hamiltonian in the " +" valley 



H, 



(+) 



efr 



(+)t 
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2 



(p(+)tp(+) + p(+)tp(+)) ,2 (p(+)t) V ,^,^p(+)ti.(+)t 



( 



+ ^1^2^2 



(+)pW 



^2^4 



.(+) 



(15) 



r 



where we introduced the kinetic momenta Pj^^ = 
F-^^'/'U- for j = 1,...,4. The Hamiltonian in the 



valley is then expressed as 



(-) 

eff 



(+)* 



eff 



The 



term yields a small correction to the hopping energy 7 
that affects the high energy bands. However, it produces 
only sub- leading corrections of order F^/j <C 1 to the 
low energy Hamiltonian (15), which are thus neglected. 



In the effective Hamiltonian (15) we also included an 



on-site energy difference A between the two layers de- 
scribing the effect of an interlayer electric field. This 
term, together with the antisymmetric component of the 
deformation potential D^^"^ has the physical effect of in- 
ducing a gap in the electronic spectrum. Thus, our re- 



sults show that such a gap is in principle realisable, with- 
out any inter-layer electric field, under different variation 
of local areas for the two layers. 



The effective Hamiltonian (15) is the main result of 



this paper. It can be used as the starting point for the 
investigation of electromechanical properties of bilayer 
graphene under arbitrary elastic deformations. Without 
loss of generality, from now on we focus on one valley 
(r = +). We first review the electronic spectrum in 
the absence of mechanical deformations, and then briefly 
discuss the effect of strain configurations inducing uni- 
form fictitious gauge fields. The corresponding effects 
on the low-energy electronic band-structure of bilayer 
graphene have been recently discussed by other authors 



7 



as well!^^^ As a further application we will then use the 



Hamiltonian (15) to describe the electron-phonon cou- 



pling and its contribution to the resistivity of suspended 
bilayer graphene. 



IV. BAND STRUCTURE WITHOUT 
DEFORMATIONS 

In the absence of mechanical deformations the effective 
Hamiltonian (p^Sl) reduces to 



(+) 

eff 



( 



if 

7 V («! + «i«2)p^ 



-A/2; 

2 



(16) 



2v^v^ \p\ 



\p? ) 



2^2 ^4 



Different cases can be analysed, according to the impor- 
tance of the various interlayer hopping terms and to the 
asymmetry in the intralayer velocities Vi and Be- 
fore discussing the most general case, we briefly review 
the features stemming from the successive introduction 
of the hopping terms, in order of descending magnitude. 

i) Most results on bilayer graphene at energies not 
too close to the Dirac point can be analysed by con- 
sidering only the dominant A2 — Bl interlayer hopping 
term 7, with equal intralayer velocities v-^ = V2 = v 
8 • 10^m-s~^, thus neglecting as well as in Eq. 



(16). In the absence of interlayer electric fields (A = 0) 
this yields a parabolic spectrum with two energy bands 
= ±^^1^1^/7 = ±|pp/2m characterised by an effec- 
tive mass m = j/2v'^ 0.05 (m^ the free electron 
mass).'^^!^ The two bands touch at one Dirac point at 
\p\ = 0. Using the complex representation of the mo- 
mentum p = \p\ exp[z^] the effective Hamiltonian can be 
written as 



(+) 



eff 



2m 



(T • (cos 2^, sin 2^, 0) 



(17) 



with cr = {(j^^cTy^cr^) the pseudospin vector of Pauli 
matrices in the sublattice space. Comparing this to 
the electronic energy, we obtain the chirality condition 
a ' (cos 2^, sin 2(/), 0) = ±1, with the upper/lower sign 
corresponding to the conduction/valence band. Thus 
electrons behave as pseudospin- 1/2 massive chiral parti- 
cles, their pseudospin winding by 47r anticlockwise when 
goes from t o 27r. This corresponds to a Berry phase 
of ^ = sQ = 27r,'^^^' where s = 1/2 is the particle pseu- 
dospin and ft = Air the solid angle enclosed by the pseu- 
dospin vector while the electronic state is transported an- 
ticlockwise in a closed loop in the 2D momentum space 
around the Dirac point. 

In the presence of a perpendicular magnetic field 
this spectrum is characterised by a doubly-degenerate 
Landau-level at zero energy (per spin and per valley) .^^^ 
Finally, an interlayer electri c field (A 7^ 0) yields an en- 
ergy gap A in the spectrum.^^ 

ii) The next step in the hierarchy of approximations is 
to include the terms proportional to ^3 10^ m • s~^ in 



Eq. ( |16| ), while still neglecting and keeping = V2 = 
V. In this case, and for A = 0, one obtains two energy 
bands with dispersion 



± 



v^p 



P 

2m 



(18) 



The two bands touch at zero energy at four Dirac points, 
obtained by imposing = 0. The four points are given 
by \p\ = as well as \p\ = 2mv^ and (/) = 7r(2n + l)/3, 
with n = 0, 1, 2, highlighting the latent symmetry of the 
honeycomb lattice. Around each Dirac point the disper- 
sion is massless. It is isotropic around the central point 
at \p\ = and anisotropic for the three satellite cones, 
as illustrated in Figs.js^ andjsjD.^^ The linearisation of 
the Hamiltonian near the new Dirac points will be pre- 
sented at the next level of approximation, where is 
also included. 

The four cones meet at energy e* = jvl/^v'^ = mvl/2 c:^ 
1.6 meV so that at higher energies the spectrum is essen- 
tially parabolic. 

As function of the electron density (i.e. of the Fermi 
energy Cp), the Fermi sea (FS) changes its shape and 
topology. In each valley, we can compactly denote the ty- 
pology of a FS characterised by electron pockets and 
hole pockets by the notation [A^^, A^^]. For |ep| > e*, 
corresponding to a concentration of electrons (or holes) 
larger than {2/7T^){mv^/h)^ 2 ■ 10^° cm-^ the FS in 
each valley is connected and topologically equivalent to 
a circle (i.e. of type [1, 0] or [0, 1]). At |ep| = e* a Lifshitz 
transition (LT) occurs and the FS shows knots that de- 
velop into four disconnected electron (hole) pockets for 
< |ep| < e* (i.e. type [4,0] or [0,4]). The single- 
particle electronic density of states (DOS) vanishes lin- 
early while approaching zero doping and exhibits a peak 
at the LT, that should lead to observable features, e.g. in 
compressibility measurements or the transport properties 
as a function of the carrier density. The topology of the 
Fermi surface close to the Lifshitz transition is presented 
in Fig. [3]d. 

In the presence of a perpendicular magnetic field, each 
massless Dirac cone yields one Landau-level at zero en- 
ergy, as in monolayer graphene, leading to a four- fold de- 
generacy per valley and spin. This scenario should result 
in quantum Hall plateaux in the transverse conductivity 
at ±8e^//i for small enough magnetic fields and close to 
zero carrier density. 

It has to be pointed out that e* is a rather small energy 
scale. The physics of the LT can thus be observed only 
in extremely clean samples at low density, so that the 
smearing due to disorder does not obscure the pertinent 
features. An alternative possibility to overcome this dif- 
ficulty would be to tune the LT to higher energies. This 
could be achieved in bilayer samples by inducing signifi- 
cant strain (see Sec.|v|). Alternatively, ABC-stacked tri- 
layer graphene shows a LT around 10 meV."^ The larger 
energy associated with this LT allows for its experimen- 
tal observation in samples with high mobility, as reported 
recently.^^ However, as the LT in trilayers originates from 
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FIG. 3. Electronic band structure of bilayer graphene with- 
out deformations. All energy scales e are in meV. a) The 
low-energy spectrum neglecting terms in V4. Four massless 
cones touch at zero energy, b) Equipotential lines for the 
panel a). Dashed lines correspond to e < e*, yielding four 
disconnected electron pockets. The thick line corresponds to 
e = e* where the LT occurs. All other continuous lines are at 
e > e*, yielding a single connected electron pocket. Dark ar- 
eas correspond to states close to zero energy while light ones 
are for higher energies, c) Band structure including V4. The 
central Dirac cones touch at zero energy, while the other three 
touch at e. The boxes [N^^N^] in different energy windows 
indicate that the FS is made out of electron pockets and 
hole pockets, d) Same as c), but with asymmetric in- 
tralayer velocities, corresponding to t2 = ti/4 = leV. This 
large asymmetry is used to stress the formation of the mini- 
gaps A. e) Magnification of the low energy spectrum in c), 
but with a small interlayer gap A < e. f) Same as in e), but 
with a larger interlayer gap A > e. 



the vertical hopping between the first and the third layer, 
it should be only weakly sensitive to deformations and it 
could not be easily tuned as in bilayers. 

iii) If we still assume v-^ = V2 = but we do not 
neglect the Al — A2 and Bl — B2 interlayer hoppings, 



taking into account v^^lA- 10 m • s ^ in Eq. (16), we 



still obtain four massless Dirac cones at low energy, but 
the two bands do not touch at the same energy.^^ Indeed, 
for A = 0, they are given by 



V 2m 



± 



P 

2m 



The four Dirac points are found at \p\ = 0, where 
the bands touch at zero energy, as well as |p| = 
2mv^/{l + vj/v'^) with (j) = 7r(2n + l)/3, where the 
bands touch at energy e = 2'yvvlv^/ {v'^ -\- vl)'^ = 

Amvlv^/ V (1 -\- vl / v'^)'^ . Since v^/v <C 1 we have 

e ^e*v^/v ~ 0.2 meV. This dispersion is illustrated 
in Fig. 

We can explicitly expand the Hamiltonian around the 
different Dirac points (described by the complex mo- 
menta pj^) by considering p = + 5p. Around = 
the linearised Hamiltonian is given by 



(+) 



v^5p 
V2,5p^ 



- 'U3 I Jp| cr • (cos (/), — sin (/), 0) 

(20) 

and describes massless chiral fermions with Berry phase 
— TT, due to the clockwise winding of the spinor for an 
anticlockwise loop of 5p. In a similar way, the expan- 
sion around one of the other Dirac points, e.g. = 
—2mv^/{l + v'l/v'^), yields 



H 



(+) 



eff 



el 





5p^ + 3z(5p^ 



(19) 



= el — \Sp\ a ■ (cos 3 sin (/), 0) , 

describing massless Dirac fermions with Berry phase tt 
and elliptical equipotential lines (see Fig. (sj)). 

As a function of the carrier density, the FS develops 
interesting structures. Due to the fact that the massless 
Dirac cones touch at different energies, the DOS never 
vanishes. For < ep < e the FS is of type [1,3]. On the 
other hand, the spectrum at energies larger than e and 
the LT remain essentially unaffected (see Fig. [s]:)). The 
critical energy for the occurrence of the LT is slightly 
renormalised to e* (1 + 2v^/v) for electron doping 

and to el c:^ — e* (1 — 2v^/v) for hole doping. 

In the presence of an external magnetic field. Landau- 
levels corresponding to massless Dirac fermions are gen- 
erated in the low energy sector. However, due to the 
energy offset e, there is only one Landau level at zero en- 
ergy (per valley and per spin) stemming from the central 
Dirac cone. This mechanism would result in plateaux of 
the Hall conductivity at ±e'^/h close to zero doping at 
low magnetic fields in extremely high mobility samples. 
However, due to the smallness of the energy scale e, the 
shift of the three satellite Dirac cones with respect to the 
central one is not observable with the present quality of 
bilayer samples. 

It is interesting to notice that a finite value of interlayer 
electric field corresponding to < A < e does not induce 
a global gap in the spectrum, but rather opens a gap in 
each individual Dirac cone. A global gap opens up only 
for A > e. Thus, in the presence of A 7^ 0, one can 
identify seven regions where the structure of the FS is as 
follows: 

1) [1, 0] for ep > e74 + V(A/2)2 + (4)2, 

2) [4, 0] for ee[i+ A/2, e/i + ,/{A/WH^i 

3) [1,0] force [max[e- A/2,A/2],e + A/2], 
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4) [1, 3] for e e [mm[i - A/2, A/2], max[? - A/2, A/2]] if 
A<i, and [0, 0] if A > i, 

5) [0,3] for e G [-A/2, min[e - A/2, A/ 2]], 

6) [0,4] for e G [e/4 - ^(A/2)2 + (e;)2, - A/2] and 

7) [0, 1] for ep < e74 - V(A/2)^ + (e* )2. 

These regions are highhghted in Fig. [Sj), and [Sf. 
iv) Finahy, we can analyse the general case in which 



and are present in Eq. (16), with Vi 7^ '^2- This 



can be of relevance for bilayer samples on a substrate, as 
the latter may induce an asymmetry in the two intralayer 
hopping energies and In this case, for A = 0, the 
low energy spectrum is given by 



^4 K 



\p\' 



(21) 



1^2 
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— 



as illustrated in Fig. [3]i. The two bands touch only at 
p = 0, where a single massless Dirac point survives while 
the other three disappear. In the experimentally relevant 
regime Vi — V2 <C ^/v]v^ ^ v the spectrum shows the 
opening of a minigap A ^ ~ \ each of the 

three satellite Dirac points. Notice that no interlayer 
term A is required in order to open these minigaps. As 
a consequence a new regime appears for e — A/2 < ep < 
e + A/2, characterised by a FS of type [1,0] similar to 
the regime 3) in the above case iii). 



V. BAND STRUCTURE WITH 
DEFORMATIONS 

In the presence of generic elastic deformations of the 
lattice, the induced fictitious gauge fields modify the elec- 
tronic low-energy Hamiltonian as in Eq. (15). While 
our formalism allows us to treat arbitrary distortions, 
as a first application we consider specific static lat- 
tice deformations which lead to uniform fictitious gauge 
fields, in analogy with those recently considered in the 
liter at ure.'^^l^ These gauge fields induce a shift in the 
electronic momenta that results in significant modifica- 
tions to the band structure. The massless Dirac cones at 
low energy drift with the deformations until they annihi- 
late at a critical value of strain. Increasing the strain fur- 
ther leaves a low energy spectrum made of two massless 
Dirac cones only. The modification to the band struc- 
ture changes the nature of the LT as well as its energy. 
Thus, controllable strain could be used to induce the LT 
at the Fermi level, with observable consequences in the 
electronic DOS and other physical characteristics. The 
deformations leading to uniform gauge fields are uniax- 
ial strain along an arbitrary direction, a rigid shift of 
one layer with respect to the other as well as a pure 
shear deformation. The first two types of distortions 
have been considered recentlyj ^^l^ although the quali- 
tative consequences in the spectrum are essentially the 



same for generic configurations leading to uniform gauge 
fields. 

A uniaxial in-plane strain along the direction = 
(cos ^, sin ^) is described by the vector u^^^ = h^^"^ = 
/i(^) = and u(^)(r) = P^^r^^O ^ P^rJ^, with = r-0, 

= and 6j_ = zxO. Here /^y and represent the 
values of the strain along the two principal directions. 
This distortion results in a uniform deformation poten- 
tial D^^"^ which is reabsorbed in a global shift of the zero 
energy, while D''^^ = 0. The corresponding gauge fields 

are given by F^^l^ = ?>a/A{dtJda){l3^^ - exp[-i2(9], 
= 0, ^3^+^ = 3aV4c(S73/Sc)(/3|| - exp[-i2l9] and 



^4^^^ = 3aV4c(S74/Sc)(/3|| - /3^)exp[-i2^]. 

Similarly, a pure in-plane shear of amplitude can 
be described by the deformation vector u^^^ = h^^^ = 
= and u^'^^(r) = C^y^x? resulting in the gauge 
3a/4(5t^/5a)(-Oexp[-z2l9], = 0, 



fields Fj 



(+) 



1=1,2 



3a'^/4.c{dj^/dc){-iC)exp[-i20] and 



(+) 



3a^/4c{djJdc){-iC) exp[-z2i9]. 

Finally, a shift of the second layer with respect to 
the first one by the amount along 6 is described 
by the deformation vector u^*^^ = u^^^ = ^a^, h^^"^ = 
= 0, resulting in the gauge fields F^^^ ^ — ^1 ^7 = 
-2{d^/dc){iaf/c, Ft^ = 3ayc{dj,/dcM)expm 
and = 3a^/c{d-f^/dc){-i^)exp[i0]. 

Quite generally, these different deformations translate 
into complex gauge fields in the Hamiltonian, which then 
affect the electronic band structure. The evolution of the 
electronic band structure under progressive strain is illus- 
trated in Fig. [4j Here we show the effect of a uniaxial 
strain of amplitude /3 = /^y — Pj_ along ^ = 0, equivalent 
to a uniform shear of amplitude C — P along 6> = 7r/4. In 
Fig.Hwe choose = and Vai,B2 = Vai,A2 = Vbi,B2 = 
1 for illustration purposes. The electronic band structure 
in the wavevector space {k^^ky) is shown in panels al) to 
a4) at different values of f3. While increasing /3 two cones 
with chirality tt and — tt approach each other until they 
annihilate at a critical strain P^-^ . Increasing P further in- 
duces the two fused cones to produce a local minimum at 
finite energy, until a second critical strain is reached. 
For P > P^2 ^tie local minimum disappears, leaving two 
massless Dirac cones at low energy. In Fig|4j panel al) 
illustrates the band structure in the regime < p < p^i 
and bl) the corresponding equipotential lines for elec- 
tronic states at positive energy. Two LT are visible at 
two different energies e^^^ > e* and €^^2 < The LT at 
e^^ separates a FS of type [1,0] for e > e^^ from a FS of 
type [3,0] for ^1^2 ^ ^ ^ Similarly, for < e < e^2 
the FS is of type [4,0]. Analogous results are obtained 
for hole doping at negative energies. 
Panels a2) and b2) present the scenario for P = P^^. Two 
Dirac cones fuse at zero energy and only one LT is left at 
^Li- Notice that the value of ej^^^ grows while increasing 
the amount of strain. The LT separates two FS of type 
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[3,0] and [1,0] for < e < e^^ and e > e^^, respectively. 
Panels a3) and b3) illustrate the regime < f3 < {3^2 
where a local minimum at finite energy e^^ survives. 
The LT separates two FS of type [3,0] and [1,0] for 
< e < and e > e^^^ respectively. A new regime 




FIG. 4. Electronic band structure in the wavevector space 
and equipotential lines for different values of strain /3 along 
^ = 0, see text. Energy (e) is expressed in meV and the 
plots are taken for = t2 = 2.47 eV and ?7ai,b2 — Vai,A2 — 
^31,32 — 1- For these parameters we get (3^1 — 2.3 • 10~^ and 
/3c2 ^ 2 • 10"^ al) Band structure for /3 = 1.5 • 10"^ bl) 
Equipotential lines for al), showing two LT at and 
(thick lines). For e > Ci^i the FS is of type [1,0] for this and 
all other panels. Dashed lines show the [3,0] FS at an energy 
between the two LT, while the thin lines exemplify a [4,0] 
FS at e < eL2- cl2) Band structure with critical strain /^d. 
b2) Equipotential lines for a2). One LT occurs at Cli, below 
which the Fermi surface is of type [3,0]. a3) Band structure 
for /3 = 4 • 10~^. The corresponding equipotential lines are 
shown in panel b3). The dashed line shows the [3,0] Fermi 
surface at < e < Cl^, while thin lines show the [2,0] FS at 
< e < Cm- a4) Band structure for /3 = 3 • 10"^. The local 
minimum disappears. The corresponding equipotential lines 
are shown in panel b4). One LT occurs at eL^. The dashed 
line shows a [2,0] FS for < e < Cl^. 



with FS of type [2,0] appears for < e < e^^. 
Finally, for {3 > {3^2 ^ single LT occurs separating FS of 
types [2,0] and [1,0] for e < e^^ and e > e^^^ respectively. 
This is illustrated in panels a4) and b4). It has to be 
noticed that in this regime of strain, the value of e^^ is 
significantly larger than e*. The ability to tune the en- 
ergy of the LT allows one to explore it at different levels 
of doping and also to partially overcome the problems in 
resolution due to disorder. 

The qualitative picture above is reproduced, essentially 
unaffected, once the strain is applied at different angles 
The precise values of the critical strains and 
13^2 depend on and on the values of the parameters 
9tsi,s'u/9^si,s'U - -Vshs'v'^si^s'i'/^si^s'V Since only the 
intralayer derivatives are known {rj^-^ = 77^2 b2 — 
it is not possible to give a quantitative estimate of the 
critical strains. However, as shown in Fig. [41 for 7?^i^^2 — 
Vai A2 — Vbi B2 — ^ typical order ofmagnitude for 
them is around 0.2 % — 2%. This is easily achieved in 
realistic suspended graphene samples. 

These parameters seem to suggest that the low energy 
spectrum in conventional suspended bilayer samples with 
high mobility is probably characterised by two massless 
cones per valley, instead of four. As a consequence, in the 
presence of an external magnetic field, the degeneracy of 
the LL at zero energy is expected to be eight, due to 
two massless cones and two spins in each valley.^ This 
picture is compatible with the recently observed integer 
quantum Hall effect at filling factors u = ±4 in bilayer 
graphene at low magnetic fields. A detailed analysis 
of uniaxially strained bilayer graphene in the presence of 
magnetic fields has been recently discussed in Ref. |17l 

As highlighted in Fig. [4j panels bl) to b4), the struc- 
ture of the Fermi surface at a given density is thus af- 
fected by strain, due to the fusion of the Dirac cones. 
As a consequence, the nature of the Lifshitz transi- 
tions is sensitive to mechanical deformations. A con- 
trollable increase of the amount of strain in the bilayer 
graphene membrane could drive the Lifshitz transition 
at the Fermi level with observable consequences on the 
electronic DOS, as illustrated in Fig. [5| This effect could 
be directly detected in compressibility measurements as 
a function of density while keeping the strain constant 
or at a fixed density while tuning the strain. The modu- 
lation of the DOS at the Fermi level could also produce 
signatures in the linear conductivity in the absence of a 
magnetic field as long as the relevant diffusion coefficient 
is smooth across the LT. 

Similar consequences to the strain-induced scenario de- 
scribed above have been found in terms of an interaction- 
induced spontaneous symmetry breaking leading to a ne- 
matic phase charac terised by two massless Dirac cones 
at low-energy !^^ * ^^l The latter scenario has been dis- 
cussed in a very recent measurement on suspended bi- 
layer graphene with extremely high mobility. 
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FIG. 5. The electronic DOS and Fermi surface for different 
values of strain jS along ^ = 0. Here we choose ?7ai,b2 — 1 ^ind 
74 = 0. a) Electronic DOS as a function of energy for /3 = 
(thick line), /S = 10~^ (dashed line), /3 = 2-10~^ (dotted line), 
^ = 5 • 10~^ (dot-dashed hne). The peaks in the DOS at the 
various LT are clearly visible, as well as the linear dependence 
on energy in the low-energy regime due to the massless Dirac 
cones. The dot-dashed line shows a step-like feature at 
associated to the local parabolic minimum in the dispersion 
(see Fig. [4] panel a3). b) Fermi surface at electron doping 
corresponding to Cp = 0.8 meV < e*. The thin line is for 
/3 = 0, yielding a [4,0] FS. The thick line shows the FS at the 
LT (ep = 61^2) for /3 1.2 • 10~^, while the dashed line shows 
the [3,0] FS at /3 :^ 3 • 10~^. c) Fermi surface at electron 
doping corresponding to Cp = 2.2 meV > e*. The thin, thick 
and dashed lines correspond to /3 = 0, /3 1.6 • 10~^ and 
P ^ 3 ■ 10~^, respectively. The LT here occurs at Cp = Cli. 



VI. PHONONIC CONTRIBUTION TO THE 
RESISTIVITY VIA FICTITIOUS GAUGE FIELD 
COUPLING 

As a further application of the gauge fields we deduced 
in Sec. |III[ we discuss the consequences of the correspond- 
ing electron-phonon coupling on the resistivity of sus- 
pended bilayer graphene. The interesting aspect to be 
pointed out is the appearance of a linear coupling be- 
tween electrons and symmetric flexural phonons {h^^"^) 

in the gauge field term F^^^^ This is in contrast to the 
case of monolayer graphene, where fiexural phonons have 
a quadratic coupling protected by symmetry with respect 
to the plane of the membrane. The contribution to the 
resistivity due to flexural modes is to be compared with 



the corresponding one due to in-plane phonons. 

In monolayers the competition is driven by a combina- 
tion of dispersion and coupling of the various phononic 
modes. In-plane ones are hard to excite (they have a 
linear dispersion) but have a strong linear coupling to 
electrons. In parallel, flexural deformations are soft (they 
have a quadratic dispersion in the absence of tension) but 
have a weak quadratic coupling. If the tension is weak, 
it has been shown that flexural phonons dominat e the 
in-plane modes as far as the resistivity is concerned. '^^'^^ 
Recent measurements of the contribution to the resistiv- 
ity due to electron-phonon scattering in suspended mono- 
layers have indeed shown the dominant contribution by 
flexural phonons. 

In bilayer graphene, the electron-phonon contribution 
to the resistivity has been considered recent ly in di- 
rect analogy with the monolayer case.l^^ESEII In-plane 
phonons as well as flexural modes with a quadratic 
coupling stemming from intralayer deformations {D^^\ 

f[^^ and ^2^^^) preserve the qualitative temperature- 
dependence of monolayer samples. However, the appear- 
ance of a linear coupling for symmetric flexural modes 
offers the opportunity of investigating soft modes (for 
weak tension) with a strong linear coupling which were 
not considered so far. Here we discuss this new contri- 
bution and analyse its consequences for transport. 



A. Phonon dispersion in bilayer graphene 

In order to study the electron-phonon contribution 
to the resistivity, we need to discuss the dispersion of 
phonons in bilayer graphene. This is easily obtained con- 
sidering the elastic Lagrangian densit};^^ 

£ = ^ /:^o^pi with 

^=1,2 

-i(K(vV))2 + r(v/i«)2) 
/:coup. = -\p.^r. (u(^) - u(2))' - \p,nl (/.(I) - /.(^))' 

where is the mass density in each layer, /i and A are the 
Lame coefficients for in-plane stretching, n is the bend- 
ing energy and F is a sample specific coefficient describ- 
ing the degree of tension induced in the membrane.!^ 
The term /^coupi i^odels a harmonic confinement for the 
sliding of one layer with respect to the other as well 
as for the modification of the interlayer distance. Re- 
cent first principle calculation^^ produced the estimates 
Q..^ ^ b - 10^2 Hz and Vt^ ^ - lO^^ jjz. The Euler- 
Lagrange equations for the elastic Lagrangian at har- 
monic level are solved in terms of the symmetric and 
antisymmetric deformations u'^^^(r) = XlqUq^^ exp[iq-r] 
and h^^^ (r) = h^^^ exp[zq • r] , with v = S^A and Uq^^ , 
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/iq^^ their Fourier transforms in the wavevector space. 
The in-plane phononic eigenmodes are given by longitu- 



dinal and transverse components i^q^'^^ = Uq"^^ • q and 

u"^'^^ = Uq^^ • q^, with q = q/|q| and q^ = z x q. The 
dispersions of in-plane as well as flexural (F) deforma- 
tions are 



1/2 



Po 



1/2 



(22) 



.1/2 



These results hold at harmonic level, while anharmonic 
corrections due to the coupling between the two layers 
would stiffen the bending coefficient k,. In realistic bi- 
layer graphene membranes symmetric flexural phonons 
would thus disperse as in the equation above, with a 
slightly renormalised bending energy. This is in anal- 
ogy with the dispersion of bending modes in carbon- 
nanotubes.ES The antisymmetric modes are gapped due 
to the harmonic interlayer couplings in C^^^^^ and do 
not give a relevant contribution to the resistivity in the 
linear transport regime. In the following we will thus 
concentrate on the symmetric flexural deformations with 
a linear coupling induced by F^'^^ and discuss their con- 
tribution in comparison with that of in-plane modes and 
of flexural modes with quadratic intralayer coupling dis- 



qu 

cussed elsewhere) ^^ * ^^ * ^-^ * We will denote the dispersion 
of symmetric flexural modes as uj^ = uj^''^^ ~ aq'^ 
(where a = ^J 4.6 • 10~^ m'^/s) for q ^ q^ and 
cjq ~ (^q^q for q <^ q^^ where q^ = \/tJk, is a sample- 
specific wavevector related to the degree of external ten- 
sion. Even in the absence of tension, anharmonic elastic 
corrections have been shown to induce a modification of 
the dispersion of flexural phonons at low energy.l^Ell 



B. The contribution of flexural phonons to the 
resistivity 

As far as the electron-phonon coupling is concerned, 
we consider the gauge field term coupling proportional to 

(t) 

F3 and focus on the high electron-density regime where 
the Fermi wave vector kp is larger than the inverse mean 
free path due to disorder and electron-phonon scattering. 
In this regime a quasiclassical Boltzmann approach to 
transport can be employed. Except for extremely clean 
bilayer samples, this condition is fulfilled for ep > e*, 
above the Lifshitz transition, where the effective elec- 
tronic dispersion is parabolic. As a consequence, the rel- 
evant electronic and electron-phonon coupling Hamilto- 
nians in one valley are given by 



(+) 



^ / pt2 

2m V 



(+) 

el— ph 



,(+)t 







An electronic eigenstate with wavevector k and en- 
ergy ej^ = Ti^k'^ /2m is described by the spinor |k) = 
(1, exp[i2^j^]), with (j)^ the angle of k with respect 
to the X axis. 

In our case the dominant coupling due to symmetric 
flexural deformations is given by ^3^^ — [dyh^^^ — 
id^h^^^)^ with ^3 = 3ac/2c (973/90) the coupling 
strength. In the Fourier space this corresponds to an 
electron-phonon coupling matrix 







-le 




(23) 



in the Dirac description, where Q is the phonon wavevec- 
tor, <I> its angle with respect to the x axis and = 
1/2 

{Ti/2MijjQj the oscillator length (M the total oscilla- 
tor mass per unit area). 

In order to calculate the resistivity, a systematic 
derivation of the Boltzmann transport equation for 
electron-phonon coupling in graphene has already been 
presented in a previous work.^^ The longitudinal resistiv- 
ity is expressed as 



m 



(24) 



where n = /cp/vr is the electronic density and 



duf 



n 



x|(k + Q|«;Q|k)|'5(ek+Q-ek) (25) 



is the scattering rate at the Fermi level in the quasi- 
elastic approximation, due to both phonon absorption 
and emission processes. Here 6 is the scattering angle 
between the electronic wave vectors k and k + Q and 
Uq = l/(exp[fico'Q//cgT] — 1) is the equilibrium Bose dis- 
tribution. The derivative of the Bose distribution im- 
plies that the relevant phonons to be considered have 
energies up to hujQ ^ k^T and their wavenumbers are 
restricted to Q ^ g^, with hujq = k^T. In this regime 
—ujQdriQ/dujQ k^T/hujQ. Implementing the on-shell 
condition due to the quasielastic approximation we ob- 
tain 

|(k + Q|^Q|k)|'=^32Q2^^sin2(3$), 
1 - cos 6> = 2 cos^ <l> 



as well as 



2m 



(5($ - $o) 



(26) 



where $o ^^o angles fulfilling the condition 



cos <I>n 



-Q/2k. As a result, rescaling Q by 2/cp, we 
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obtain 

1 Wmngl 



(27) 



^0 



/ ,2 



Different regimes appear as a function of temperature 
and of tension. The low temperature regime qrp/2kp <C 1 
corresponds to T <C ^^^^ ^bg ~ ^^2k /^b 

Bloch-Griineisen temperature, and is dominated by small 
angle scattering which give little contribution to the re- 
sistivity due to the term 1 — cos6>. In addition, for weak 
tension, flexural phonons in bilayer graphene are charac- 
terised by TgQ 0.4 nK, with h the density expressed in 
units of 10^^ cm~^. As a consequence, for typical param- 
eters of relevance to experiments, electron-phonon scat- 
tering yields a significant contribution to the resistivity 
only in the high-temperature regime, T ^ ^bG' where 
Qrp ^ 2/cp. In this regime, for weak tension {q^ <C 2/cp) 
the resistivity takes the value 



r weak— tens 



h rn^gl k^T 
^ ' 87r2?iVn2 



(28) 



10" 



T[K] 



In the opposite regime {q^ 2kp) dominated by tension 
we obtain the sample specific result 



r tens 



2A^p 



2 ' Pweak— tens Pweak— tens * (^^) 



As in the case of monolayer graphene, external tension 
stiffens the flexural phonons and reduces their density of 
states without affecting their coupling, thereby suppress- 
ing the contribution to the resistivity. 

For completeness, in the low temperature regime T <C 
TgQ, we have the two results (for weak tension q^ <C 
and strong tension q^ ^ qT) 



„(^«Tbg) 
^weak— tens 



r tens 



h 



Snp^h'^ {haq^{2kp) f 



(30) 
(31) 



As far as flexural phonons are concerned, the contribu- 
tion to the resistivity we just considered has to be com- 
pared with the corresponding one stemming from the 
usual quadratic coupling (induced via the deformation 
potential D^^'> as well as the intra-layer gauge fields F^'^^ 

and i^2^^)- absence of tension, and for T ^bG' 

this has been estimated^^ to be 



m 



(T»Tbg) _ _ 



(32) 



with the screened deformation potential coupling con- 
stant. As a consequence, the ratio of the two contribu- 
tions is 



(^»^bg) 
'weak— tens 

Pf 



5-W 



9l 



500 



glT[K] T[K] 



(33) 



where the last estimate is obtained with the approxima- 
tions ^3 —73/2 — 0.15 eV (corresponding to rjj^^ ^2 — 
1) and 3.5 eV. As a result, in the absence of tension 
and up to room temperature the linear coupling for flexu- 
ral phonons mediated by the interlayer gauge fields yields 
a dominant contribution to the resistivity with respect to 
the quadratic one. In this case, the temperature depen- 
dent resistivity in suspended bilayer samples is expected 
to show a linear-T dependence even in the absence of ten- 
sion, in contrast to the monolayer case. Tension would 
then suppress the contribution due to flexural phonons 
in favour of the in-plane ones. The latter have been re- 
cently discussed in bilayer graphene.'^^'^^'^The resulting 
contribution to the resistivity is 



h 



10-^ nT[K\c^ 2.6- 10-^ hT[K\n (34) 



showing a linear-T dependence as for p^^^^ens 

Ptens • The contribution due to in-plane phonons at 
room temperature is thus supposed to be about 101^ 
at n = 1, significantly less than flexural modes with 
weak tension. The critical tension F needed to suppress 

PtenT^^^^ with respect to corresponds to an in- 

duced strain of order 3-10~^/n^. Thus, as in monolayers, 
only samples with a small intrinsic strain have a chance 
to show signatures of flexural phonons. Besides their 
magnitude, the in-plane and flexural contributions show 
the same temperature dependence, but behave differently 
with respect to density. Flexural modes would indeed 
yield a contribution to the resistivity which depends on 
the density as for weak tension and as in pres- 
ence of strong tension. In contrast, in-plane modes result 
in a linear dependence of the resistivity on n. If exper- 
iments show a linear dependence of the electron-phonon 
resistivity on T, the corresponding density dependence 
would discriminate which phononic branch is dominant, 
and if tension is of special relevance. 



VII. CONCLUSIONS 

In this paper we analysed the consequences of generic 
elastic deformations on the electronic properties of bi- 
layer graphene membranes. We deduced the fictitious 
gauge fields, induced by arbitrary distortions, in the elec- 
tronic Dirac Hamiltonian and focused on the low-energy 
effective theory describing the two quasi-degenerate elec- 
tronic bands close to zero energy. 

As a first application we analysed the effect of static 
deformations yielding uniform fictitious gauge fields on 
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the four massless Dirac cones induced by trigonal warp- 
ing at low energy. An increasing degree of strain leads to 
the annihilation of two Dirac points, leaving the two re- 
maining cones as the relevant low energy band-structure. 
This effect has direct consequences on the quantisation 
of the Hall conductivity close to zero density.!^ In par- 
allel, strain allows one to tune the topological Lifshitz 
transition at the Fermi level, with additional observable 
signatures in the single particle density of states. 

We further considered the derived gauge fields as an 
electron-phonon coupling mechanism and discussed the 
consequent contribution to the temperature-dependent 
resistivity. We pointed out the appearance of a lin- 
ear coupling between electrons and symmetric fiexural 
phonons due to the inequivalence of the two layers and 
the non- vertical interlayer hopping processes. For sus- 
pended bilayer membranes with low tension, this results 



in a phononic contribution to the resistivity larger than 
the one stemming from the conventional quadratic cou- 
pling for fiexural modes as well as that due to in-plane 
phonons. 

Our investigation sets the basis for future works on 
the electromechanical properties of suspended bilayer 
graphene membranes, including applications in strain- 
engineering and in the creation of fictitious magnetic 
fields. 
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